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The  US  Air  Force  Environmental  Technical  Applications  Center  (USAFETAC)  pre¬ 
pared  this  global  atlas  of  mean  sea-level  pressure  and  its  standard  deviation  for 
Headquarters,  Strategic  Air  Command  (SAC)  under  USAFETAC  Project  2242.  It  repre¬ 
sents  a  first  for  USAFETAC  and  indeed  for  Air  Weather  Service  (AWS)  in  that  it  is 
the  initial  application  of  a  new  computer-graphics  capability. 

The  maps  are  based  on  the  UAPIP  data  set  which  is  derived  from  the  summarized 
Air  Force  Global  Weather  Central  (AFGWC)  upper-air  analyses.  Although  primarily 
an  upper-air  data  set,  the  UAPIP  was  used  because  it  was  the  only  gridded  data 
set  available  at  the  time  offering  mean  sea-level  pressure  data  over  the  entire 
globe . 

The  plotting  software  used  to  produce  the  maps  is  a  modified  version  of  the 
OSUCON  package  developed  at  Oklahoma  State  University.  It  is  discussed  below. 

In  this  atlas,  the  global  patterns  of  mean  sea-level  pressure  and  the  stand¬ 
ard  deviation  of  the  pressure  are  presented  by  month  and  hemispheric  quadrant. 
Each  pressure  chart  and  its  corresponding  standard  deviation  chart  are  found  on 
facing  pages. 


2.  The  UAPIP  Data  Set 

The  UAPIP  is  a  gridded  global  data  set  containing  pressure,  D-value,  tempera¬ 
ture,  dew  point,  wind,  and  density  information.  It  is  derived  from  the  Summar¬ 
ized  Analysis  Data  Set  (SADS)  which  is  built  from  the  AFGWC  upper-air  analyses. 
For  each  hemisphere  (North  and  South)  the  UAPIP  is  presented  on  a  square  65x65 
point  grid  ( see  Figures  1  and  2 ) . 

Let  us  take  a  quick  look  at  the  development  of  the  UAPIP  data  set.  We  will 
briefly  trace  the  processes  starting  with  an  observation  arriving  at  AFGWC  and 
ending  with  the  climatological  data  formatted  into  the  UAPIP.  For  greater  detail 
the  reader  is  referred  to  the  references  in  Section  5. 

The  analysis  and  forecast  models  at  AFGWC  use  gridded  data.  Three  grids  are 
used,  the  Northern  Hemisphere  (NH)  octagon,  the  Southern  Hemisphere  (SH)  octagon, 
and  the  tropical  strip  or  grid  ( see  Hoke ,  et  al . ,  1981 ) .  There  are  approximately 
5000  reporting  stations  in  the  NH  octagon,  250  in  the  SH  octagon,  and  1150  in  the 
tropical  grid.  The  quantity  and  quality  of  the  data  from  these  roughly  6400 
meteorological  reporting  stations  world-wide  varies  with  geographical  area  and 
reporting  time. 

Radiosonde  observations  (raobs),  pilot-balloon  observations  (pibals),  air¬ 
craft  observations,  surface  observations,  satellite  derived  data,  and  bogus  data 
(see  Tarbell  and  Hoke,  1979)  are  used  in  the  upper-air  analyses.  The  number  of 
inputs  may  vary,  but  there  are  usually  700  raobs,  1100  pibals,  700  aircraft 
observations,  and  5000  surface  observations.  Satellite  derived  data  consist  of 
approximately  200  soundings  and  several  hundred  cloud  track  wind  observations. 
The  winds  are  generally  based  on  the  movement  of  low  clouds,  such  as  cumulus,  and 
high-level  clouds,  such  as  cirrus.  Bogus  data  are  input  by  the  AFGWC  Forecasting 
Services  Division  forecasters  when  they  determine  that  the  first  guess  field  is 
not  satisfactory  in  data  sparse  regions.  The  amount  of  bogus  data  input  may 
vary  from  none  to  over  100  data  points. 

AFGWC  uses  the  Hough  Global  Spectral  Analysis  Model  (HUFANL,  see  Tarbell  and 
Hoke,  1979)  to  convert  station-specific  data  to  gridded  data.  This  model  re¬ 
solves  the  global  atmosphere  with  sines  and  cosines  in  longitude.  Hough  functions 
in  latitude,  and  empirically  derived  orthogonal  functions  in  the  vertical.  This 
technique  uses  meteorological  type  (barotropic)  waves  to  resolve  the  atmosphere. 
The  smallest  resolvable  atmospheric  wave  is  15°  in  longitude.  One  might  also 
note  that  the  HUFANL  does  not  really  analyze  for  surface  pressure  and  surface 
temperature.  The  surface,  or  mean  sea  level,  values  are  actually  derived  from 
the  values  computed  for  the  1000  mb  surface.  For  terrain  above  1000  mb,  the 
1000-mb  D-values,  or  heights,  are  calculated  using  the  hypsometric  equation. 


Figure  1.  Northern  Hemisphere  65x65  Grid. 
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Once  the  AFGWC  coarse-mesh  upper-air  analyses  are  complete,  they  are  sent  to 
OL-A,  USAFETAC,  Asheville,  North  Carolina,  for  inclusion  in  the  Air  Weather  S~_v- 
ice  Climatological  Data  Base.  OL-A  then  uses  these  analyses  to  build  the  SADS. 

First,  each  meteorological  element  is  summarized  for  each  of  the  16  standard 
radiosonde  levels  by  year,  month,  and  hour  (00Z  and  12Z)  for  each  of  the  three 
analyses  (NH  octagon,  SH  octagon,  and  tropical  strip).  Table  1  shows  which  ele¬ 
ments  and  which  statistics  are  gathered  at  each  level.  Note  here,  and  while 
looking  at  Table  1,  that  the  SADS  provides  both  00Z  and  12Z  data.  The  UAPIP 
derived  from  the  SADS  is  an  "all-hours11  climatology. 

The  next  step  in  the  development  of  the  SADS  is  to  prepare  two  whole- 
hemisphere  grids  using  data  from  the  NH  octagon,  SH  octagon,  and  tropical  strip. 
This  gives  a  summarized  data  set  for  each  hemisphere.  In  combining  each  octagon 
with  the  tropical  strip  to  obtain  the  hemispheric  data  set,  the  summarized  tropi¬ 
cal  strip  data  on  a  19x72  Mercator  grid  are  translated  to  the  polar  stereographic 
projection.  This  gives  an  extension  of  the  octagon  grid  to  a  whole-hemispheric 
grid  extending  just  beyond  the  equator.  The  values  for  each  tropical-region  grid 
point  from  the  Equator  to  40.97°  on  the  polar  stereographic  projection  are  ob¬ 
tained  by  bilinear  interpolation  between  the  four  surrounding  grid  points  on  the 
mercator  projection.  The  parameter  values  for  the  octagon  grid  points  poleward 
of  40.97°  are  combined  with  those  derived  for  grid  points  from  the  tropical  strip 
to  form  a  hemispheric  polar  stereographic  65x65  grid.  The  65x65  grid  for  each 
hemisphere  thus  provides  complete  hemispheric  coverage  by  virtue  of  its  extension 
beyond  the  Equator.  This  data  set  in  the  SADS  format  is  maintained  by  year- 
month-hour  and  is  the  OL-A,  USAFETAC  summarized  upper-air  analysis  data  set. 

The  UAPIP  is  one  of  the  climatological  summaries  derived  from  the  SADS.  The 
UAPIP  climatology  contains  summarized  means  and  standard  deviations  of  sea-level 
pressure,  D-value,  temperature,  dew-point  depression,  and  the  u  and  v  wind  com¬ 
ponents.  Means  and  standard  deviations  are  calculated  according  to 


Mean  =  x  =  Ix/n  (1) 

Standard  Deviation  =  ox  =  {[nix2  -  (Zx)2]/n(n  -  l)}1*  (2) 


where  x  is  the  element  being  summarized  and  n  is  the  number  of  elements  (number 
of  observations)  summed. 

It  is  important  to  note  here  that  the  means  and  standard  deviations  in  the 
UAPIP,  and  hence  the  means  and  standard  deviations  plotted  in  the  atlas,  are  cal¬ 
culated  from  all  available  data.  This  contrasts  with  some  other  currently  avail¬ 
able  atlases  which  show  the  average  of  several  monthly  means  as  the  mean  for  that 
month  (e.g.,  taking  the  average  of  the  mean  temperature  in  January  1960,  1961, 
...,  and  1969,  and  calling  that  number  the  mean  temperature  for  January).  Where¬ 
as  the  mean  of  means  may  approximate  the  true  mean,  a  far  more  heinous  transgres¬ 
sion  is  to  calculate  a  standard  deviation  based  on  only,  say,  ten  mean  values, 
and  then  call  that  the  standard  deviation  of  the  population  when  in  fact  it  is 
simply  the  standard  deviation  of  the  monthly  means.  Any  relation  between  the  two 
will  be  purely  coincidental. 


3.  The  OSUCON  Plotting  Routine 

To  plot  the  maps  of  mean  sea  level  pressure  and  its  standard  deviation,  we 
used  the  OSUCON  software  package.  This  is  a  general  purpose  system  for  generat¬ 
ing  plotted  maps  and  charts.  At  USAFETAC,  the  OSUCON  package  runs  on  the  IBM 
4341  Computer.  A  plot  transfer  tape  is  produced  which  is  then  input  to  the 
DEC  PDP  11/45  computer  where  finished  maps  are  produced  by  a  VERSATEC  2160A 
raster  plotter. 
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Table  1.  Summarized  Analysis  Data  Set  (SADS). 


Element* 

Statistics* 

Pressure  Level  (millibars) 

Pressure  -  P  (mb) 

IP, 

IP2, 

N 

SFC  (Sea  level) 

D-Value  -  h  (m) 

Ih, 

Ih2 , 

N 

1000,  850,  700,  500, 

300,  250,  200,  150,  100, 
50,  30,  20,  10 

400, 

70, 

Temperature  -  T  (K) 

IT, 

IT2, 

N 

1000,  850,  700,  500, 

300,  250,  200,  150,  100, 
50,  30,  20,  10 

400, 

70, 

Dew-Point  Depression  -  Tj 
(K)  d 

*Td' 

IT|, 

N 

SFC,  850,  700,  500,  400, 

300 

Zonal  Wind  Component  -  U 
(m  sec" 1 ) 

IU, 

IU2, 

N 

1000,  850,  700,  500, 

300,  250,  200,  150,  100, 
50,  30,  20,  10 

400, 

70, 

Meridional  Wind 

Component  -  V 
(m  sec" 1 ) 

IV, 

IV2, 

N 

1000,  850,  700,  500, 

300,  250,  200,  150, 

70,  50,  30,  20,  10 

400, 

100, 

Cross  Product  of 

Wind  Component  -  UV 
(m  sec"1 ) 

IUV 

1000,  850,  700,  500, 

300,  250,  200,  150, 

70,  50,  30,  20,  10 

400, 

100, 

Density  -  p  (g  m"3) 

Ip. 

IP2. 

N 

1000,  850,  700,  500, 

300,  250,  200,  150,  100 
70,  50,  30,  20,  10 

400, 

♦Elements  and  associated  statistics  are  maintained  for  each  coordinate  grid  point 
for  each  period  of  record  year-month-hour  for  hours  002  and  122.  Parameter  N  is 
the  number  of  occurences  (by  element)  which  are  summarized.  The  I-J  grid  spacing 
is  approximately  200  nautical  miles  at  60°N  and  60°S. 

********** 

When  preparing  the  maps  for  this  altas,  we  specified  the  data  to  be  plotted, 
the  contour  interval  to  be  used,  and  the  area  of  and  the  orientation  of  the  map. 
The  UAPIP  data  were  read  from  disk  and  the  maps  produced  in  a  straight-forward 
manner . 

A  quick  look  through  some  of  the  maps  shows  some  of  the  peculiarities  and 
idiosyncrasies  of  OSUCON.  In  particular,  note  the  tropical  regions.  Meteorolo¬ 
gically  speaking  here  we  find  rather  flat  fields.  That  is  to  say  the  pressure 
gradients  are  small  and  pressure  change  (hence  standard  deviation)  is  similarly 
small.  Because  of  the  flatness  of  the  fields  to  be  contoured,  single  grid-point 
values  can  generate  high  and  low  pressure  centers  and  centers  of  maximum  and  min¬ 
imum  standard  deviation  that  are  not  really  meaningful. 

Another  quirk  is  the  occasional  appearance  of  "single-line”  ridges  and 
troughs  and  of  right-angle  kinks  in  the  contours.  These  arise  when  the  contour¬ 
ing  software  on  the  computer  attempts  to  analyze  very  flat  fields  and  fields 
where  little  data  is  available  as  explained  below. 

Flat  fields  give  rise  to  isolated  line  segments.  What  happens  is  that  a  few 
data  points  are  high  (or  low)  enough  to  warrant  a  closed  contour.  The  computer 
calculates  and  the  plotter  draws  an  isopleth  around  an  extremely  narrow  area.  As 
a  result,  what  is  really  a  closed  contour  appears  as  an  isolated  line  segment. 

Sparse  data  (especially  near  the  Equator  since  the  gridded  data  extend  only 
to  the  equator  and  do  not  cover  the  entire  65  x  65  grid)  gives  rise  to  sharp¬ 
angled  bends  and  kinks  on  the  maps.  In  this  case  there  is  simply  insufficient 
data  for  the  computer  to  generate  a  smooth  contour. 


4.  Map  Use  and  Comments 


The  maps  in  this  atlas  show  mean  sea-level  pressure  for  each  month  and  the 
standard  deviation  of  the  mean  sea-level  pressure  by  month.  These  maps  are  based 
on  a  6-year  period  of  record  (1973-1978). 

The  maps  show  an  averaged  pressure  pattern.  Rarely  if  ever  will  a  particular 
day's  surface  chart  look  like  the  applicable  monthly  mean  chart.  Areas  of  mean 
high  pressure  depict  regions  that  are  either  dominated  by  high  pressure  cells 
most  of  the  time  or  experience  extremely  high  pressure  values  when  transient  sys¬ 
tems  pass.  Similar  arguments  apply  to  areas  of  low  mean  monthly  pressure. 

Regions  with  relatively  high  standard  deviations  are  regions  where  the  normal 
range  of  pressure  values  is  large.  These  areas  are  dynamically  active,  experienc¬ 
ing  the  somewhat  frequent  passage  of  highs  and  lows.  Regions  of  relatively  low 
standard  deviation  reflect  the  rather  stable,  constant  weather  regimes  in  those 
areas . 

The  prime  usefulness  of  these  maps  is  to  assess  mean  sea-level  pressure  pat¬ 
terns  and  pressure  values  over  areas  and  to  estimate  pressure  values  and  ranges 
of  pressure  values  at  locations  where  a  point  summary  such  as  a  Revised  Uniform 
Summary  of  Surface  Weather  Observations  (RUSSWO)  is  not  available. 

When  using  the  maps  of  this  atlas  certain  cautions  must  be  exercised.  First, 
keep  in  mind  the  short  (especially  for  climatological  studies)  period  of  record 
used.  Next,  realize  the  different  amounts  and  varying  quality  of  the  data  that 
went  into  the  Summarized  Analysis  Data  Set  (SADS).  Finally,  do  not  try  to  say 
more  than  the  data  allow. 

Note  that  on  the  maps  north  latitude  (Northern  Hemisphere)  is  positive  and 
southern  latitude  (Southern  Hemisphere)  is  negative.  For  example,  40°N  is  40  00 
on  the  maps  and  40°S  is  -40  00 . 

Also  note  that  east  longitude  is  positive  and  west  longitude  is  negative. 
For  example,  50°w  is  -50  00  and  50°E  is  50  00. 
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